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Abstract 

We explore the Gell-Mann-Okubo mass relation among the octet baryons using fully-dynamical, 
mixed-action (domain-wall on rooted-staggered) lattice QCD calculations at a lattice spacing of 
h ~ 0.125 fm and pion masses of m^ ~ 290 MeV, 350 MeV, 490 MeV and 590 MeV. Deviations 
from the Gell-Mann-Okubo mass relation are found to be small at each quark mass. 



I. INTRODUCTION 

A decade before the formulation of Quantum Chromo-Dynamics (QCD) as the theory of the 
strong interactions, Gell-Mann [l| and Ne'eman |y] proposed that the low-lying hadrons could 
be embedded into irreducible representations of the flavor group SU{3), with a fundamental 
representation containing three quarks: up, down and strange. This provided a relatively 
simple theory for understanding the plethora of strongly-interacting particles that had been 
seen at that time, although it was a radical proposal as the quarks had not been seen in 
isolation. Now that QCD is known to be the theory of the strong interaction, classification of 
the strongly-interacting particles by their transformation properties under the flavor group 
SU{3) (and the isospin subgroup SU{2)) follows naturally. Of course, the flavor symmetries 
are only approximate, as they are broken by the differences in the quark masses. The 
impact that these mass differences have on hadronic spectra and interactions is determined 
by how they compare to the scale of strong interactions. The SU(2) symmetry of isospin 
is a very-good symmetry of nature, being violated only at the percent level by both the 
electromagnetic interactions and the mass difference between the up and down quarks. The 
approximate SU{3) symmetry is violated typically at the ~ 30% level due to the large mass 
difference between the strange quark and the up and down quarks as compared to the scale 
of the strong interaction. However, there are instances where SU{3) symmetry works much 
better than this, as in the case of matrix elements of an SU{3) charge operator between 
states in the same irreducible representation, which are protected by the AdemoUo-Gatto 
theorem J3|. 

In addition to providing a well-defined scheme for understanding the spectrum of hadrons, 
flavor symmetry also provides a vital tool for the classification of operators, such as those 
responsible for nonleptonic decays. It is well established that the octet component of 
the nonleptonic weak interactions dominates over the 27 component, a feature known as 
octet dominance. This is a generic feature of strong interactions: the contributions from 
higher- dimensional representations of SU{3) are suppressed compared to those from lower- 
dimensional representations. The origin of such suppressions is not well understood, but 
is found in perturbative calculations, and also by consideration of the large-A^c* limit of 
QCD. The Gell-Mann-Okubo mass relation among the baryons is a classic example of this 
feature of QCD. When all three quark masses are equal, the light-quark mass matrix trans- 
forms as a singlet under flavor SU{3), and as the strong interactions are also flavor singlets, 
higher-dimensional representations of SU{3) are not induced, and 

Mn = Mj^ = Ma = Ms , (1) 

where N = p,n denotes the nucleon isodoublet, S = S+, S°, S^ denotes the E isotriplet, 
and S = S^, S° denotes the H isodoublet. In the limit of exact isospin symmetry, but broken 
SU{3) symmetry, the light-quark mass matrix takes the form 

/m \ 
niq = \ m \ , (2) 

V m, ; 

where m denotes the mass of the up and down quarks, while rus denotes the mass of the 
strange quark. The mass matrix rriq transforms as an 8 © 1 under SU{3) transformations. 
A single insertion of niq in an octet baryon mass, or, more generally, an insertion in the 
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0.007 


293.1 ±1.5 


1.983 ±0.020 


2.4 ±1.8 


0.8128 ± 0.0082 


3.0 ±2.3 


0.010 


354.0 ±0.9 


2.332 ±0.012 


3.3 ±1.0 


0.8218 ± 0.0080 


4.1 ±1.3 


0.020 


492.6 ±1.1 


3.077 ±0.014 


0.34 ±0.56 


0.8724 ± 0.0055 


0.39 ±0.64 


0.030 


591.8 ±0.9 


3.588 ±0.011 


0.07 ±0.42 


0.9162 ± 0.0041 


0.07 ±0.46 



TABLE I: Violations of the Gell-Mann-Okubo mass relation computed with lattice QCD. The first 
and second columns list the sea-quark masses employed in this work and the corresponding pion 
masses in MeV, respectively. We have used h = 0.125 fm to set the scale. The fourth column 
corresponds to T defined in eq. ^, measured in lattice units (l.u.), which is also the numerator of 
eg. TO. The fifth column corresponds to the location of the centroid of the lowest-lying baryon octet 
(also in lattice units), which is the denominator of eq. TO. All errors quoted are statistical. The 
systematic errors due to fitting are all significantly smaller than the statistical errors and thus are 
not shown. 



baryon mass that transforms as an 8 under SU{3) ^, splits the masses of the baryons, but 
leaves one relation among the masses intact, 



T = 



Ma + |Me - IMn 



iMs 







(3) 



the Gell-Mann-Okubo (GMO) mass relation [l|,|^. Multiple insertions oirrig, either through 
operators in the low-energy effective field theory, or through loops involving the light mesons, 
or both, will lead to violations of this relation. As T vanishes for SU{3) breaking that 
transforms as 8 © 1 , a non-zero value of T can result only from contributions transforming 
as 27. 

The typical scale of SU{3) breaking is set by the differences in the squares of the meson 
masses, 



Q 



m 



K 



mz 



^l 



0.23 



(4) 



where A^ ~ 1 GeV. The mass differences between baryons in the octet that differ in 
strangeness by one unit are then expected to scale as AM ~ QA^, and one naively ex- 
pects violations of the GMO relation to be T ~ Q'^^x ~ ^^ MeV. Experimentally, at 
{m^/fnY''^^ = 1-045 , one finds that T'^^p* = 8.76 ± 0.08 MeV 2, where the isomultiplet 
averages of the masses have been used in eq. Q. Some theoretical understanding of why 
the GMO mass relation works better than naive expectations would suggest can be gained 
by considering the large- Ac limit of QCD, where it is found that violations to this relation 
are further suppressed by a factor of 1/Nc l^ giving T ~ Q'^A^/Nc ~ 17 MeV, which is 
close to the experimental value. 

Given the smallness of the violation of the GMO relation found experimentally, it is 
natural to ask whether this smallness persists at other values of the light-quark masses. 



^ Such a contribution will arise from n insertions of niq, i.e. resulting from the 8's in the product [1 © 8] (g) 

[1©8](8)...(8) [1®8]. 
^ The electromagnetic contribution to the violation of the GMO mass relation is expected to be an order 

of magnitude less than T°^p^. 
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or is no more than an accident for the physical values of the quark masses. In this work 
we address this issue by computing deviations from the GMO mass relation at a num- 
ber of different light-quark masses. Following the LHP collaboration (LHPC) ja |7|, our 
computation is a hybrid lattice QCD calculation using domain-wall valence quarks from a 
smeared-source on three sets oi Nf = 2 + 1 asqtad-improved [8|, |9| MILC configurations 
generated with rooted-staggered sea quarks [lO|. In the generation of the MILC configu- 
rations, the strange-quark mass was fixed near its physical value, hmg = 0.050, (where b 
is the lattice spacing) determined by the mass of hadrons containing strange quarks. The 
two light quarks in the four sets of configurations are degenerate (isospin-symmetric), with 
masses bmi = 0.007, 0.010, 0.020 and 0.030. Some of the domain-wall valence propagators 
were previously generated by LHPC on each of these sets of lattices. The domain-wall height 
is M5 = 1.7 and the extent of the extra dimension is L5 = 16. The parameters used to gen- 
erate the light-quark propagators have been "tuned" so that the mass of the pion computed 
with the domain-wall propagators is equal (to few-percent precision) to that of the lightest 
staggered pion computed with the same parameters as the gauge configurations [lO|. The 
MILC lattices were HYP-blocked [llj and Dirichlet boundary conditions were used to reduce 
the time extent of the MILC lattices from 64 to 32 time-slices in order to save time in prop- 
agator generation. Various parts of the lattice were employed to generate multiple sets of 
propagators on each lattice. We analyzed three sets of correlation functions on 564 lattices 
with bnii = 0.030, 486 lattices with bmi = 0.020, 658 lattices with bmi = 0.010, and 468 
lattices with brrii = 0.007. The lattice calculations were performed with the Chroma soft- 
ware sui 
(JLab). 

We denote as Cx{t) the correlator as a function of time slice for baryon X, with large-time 
behavior 

Cxit) ^ e-^-* . (5) 

We found that the deviations from the GMO mass relation were determined most precisely 
by forming the products and ratios of correlation functions associated with each of the 
baryons in the octet, 

^GMO/^N _ CAJt) CY,{t) -(MA+ME/3-2Miv/3-2MH/3)t /r\ 

The resulting correlation functions are shown in fig. ^and corresponding effective mass plots 
are shown in fig. El It is convenient to define a normalization factor with which to measure 
deviations from the GMO mass relation. We normalize the deviation to the centroid of the 
baryon octet, and define the fractional violation of the GMO mass relation, (5gmo? to be 

_ Ma + |Ms - iMj, - IM^ 
^^^° - |Ma + |M, + \Mx + IMh • ^^^ 

The experimental deviation translates into: 5[^mo ~ 0.00761 ± 0.00007. The results of our 
fully-dynamical lattice QCD calculation are shown in Table Ql and are plotted in fig. |21 The 
computation shows that the deviations from the GMO mass relation are small, and in fact 
less than the experimental value, at all quark masses accessible to us. It is now quite clear 
that the small violation of the GMO relation observed experimentally is not an accident 
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FIG. 1: The logarithm of the correlation function G {t) for the four quark masses employed in 
this work. Each has been off-set vertically for display purposes. 



due to the physical values of the light-quark masses, but is a generic feature of the baryon 
spectrum. 

Heavy-baryon chiral perturbation theory IJ, |l5| (HBxPT) can be used to compute cor- 
rections to the GMO mass relation. The leading-order correction arises at one-loop in the 
chiral expansion, and is dictated by the mass differences between the pseudo-Goldstone 
bosons, the vr, K and i]. It is found that |l6l.ll7lll8l 



Ma + iMs - IMn - |Mh 



247r/2 



_p2 _ 2J^2^ (^ 4^3^ _ 3^3 



ni" 



-— ( 4Fk - 3F, 



F 



where the function Fc = F{mc, A, /x) is 

F{m,A,fi) = (m^-A^) [ VA^ - m^ log 
1 



'A- VA2-m2 + 



le 



A + VA2 - m2 + 



i€ . 



Alog^ 



-Am log 



m 



(8) 



2 \fi 

and we hold the Delta- Nucleon mass-splitting. A, fixed at its physical value, A"^"^^^ 



(9) 
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FIG. 2: Effective mass plots for the correlation function G'^^'^{t) at the four quark masses employed 
in this work. 
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FIG. 3: Deviations from the Gell-Mann-Okubo mass relation for the lowest-lying baryon octet. The 
red diamond corresponds to the experimental value. The brown band corresponds to the leading- 
order one-loop SU{3) HBxPT prediction, and its expected ~ 30% uncertainty. The black data 
points are the results of our fully- dynamical lattice QCD calculation at pion masses of tti^ ~ 
293 MeV, 354 MeV, 493 MeV and 592 MeV. 



^. To the order we are working, the fi independence of the baryon GMO re- 
(jH)). is protected by the meson GMO relation, and thus no counterterms are 



nrrip. 



293 MeV 
lation, eq. 

required. For vanishing decuplet-octet mass sphtting, the function F tends to Fp 
The constants V and J-" are the axial couplings between the octet baryons and the octet 
of pseudo-Goldstone bosons, while the constant C is the axial coupling among the decu- 
plet of baryon resonances, the octet of baryons and the octet of pseudo-Goldstone bosons. 
At tree-level these couplings take the values [19j JF = 0.5, V = 0.8 and C = 1.5, while 
at one-loop level they take the values [19| J-" = 0.4, V = 0.6 and C = 1.2. In fig. 01 we 
have shown (the brown band) the one-loop prediction of HBxPT, which is seen to agree, 
within uncertainties, with both the experimental value, and also the results of our lattice 
calculation. In obtaining this curve, we have divided the HBxPT expression of eq. (jH)) by 
the experimental value of the octet centroid. As we have not computed the mass of the 
?7-meson on the lattice due to the problems associated with disconnected diagrams, we have 



The quark-mass dependence of the decuplet parameters is a higher-order effect. 



used the Gell-Mann-Okubo mass relation among the square of the meson masses as input 
for TJirj into the one-loop chiral expression in eq. (jHI), and used our lattice data for the tt 
and K parameters. Expressions for the octet baryon masses at higher orders in the chiral 



expansion have been computed [17|, |20|, |21| . However, with lattice calculations at only four 
pion masses, and the heaviest likely outside the range of validity of the chiral theory, there 
is little point in attempting to fit such expressions. 

As with any lattice QCD calculation, there are additional systematic uncertainties that 
must be estimated or removed. The finite-volume effects on the baryon masses are estimated 
to be negligible. The MILC lattices used in this work have spatial extent L ~ 2.5 fm and 
finite- volume chiral perturbation theory for baryons J22] indicates that for the meson masses 
used, the finite-volume effects are very small. The source of uncertainty in this calculation 
that we have the least control over is that due to the finite lattice spacing, b, however such 
corrections are expected to be small. The lattice spacing transforms as a singlet under flavor 
transformations, and consequently, it does not contribute at leading order to violations of the 
GMO relation. Even corrections of the form 0{mqh'^) will not contribute as they transform 
as an 8© 1. However, we expect contributions of the form f Am\- — 3m^ — m^ j 6^, which is 
an element of the 27. Further, there will also be exponentially-suppressed 0{h) contributions 
multiplying objects that transform as a 27. The centroids of the baryon octet will experience 
0(6^) shifts, as well as the exponentially suppressed 0{h) contributions. A full analysis of the 
chiral corrections at finite lattice spacing should be performed in the mixed-action theory [23| 



when resources to calculate at another lattice spacing materialize. While there have been 



other quenched QCD [2J, |25|, |26|, |27| and fuU-QCD i28„ ^ ^, ^J studies of octet baryon 
masses, to our knowledge there have been no other fuU-QCD studies of the complete baryon 
octet in the chiral regime, prior to this work. 

II. CONCLUSIONS 

In this paper we have explored violations of the Gell-Mann-Okubo mass relation among 
baryons in the lowest-lying octet using fully-dynamical lattice QCD. With domain-wall va- 
lence quarks on rooted staggered MILC configurations, we have calculated the deviations to 
the GMO mass relation at four different pion masses. We find that the violations are con- 
sistent with the size of the experimental value, and vanish toward the limit of exact S'f/(3) 
symmetry as required. Therefore, we conclude that the smallness of the experimentally- 
observed violation of the GMO mass relation is not an accident. This is consistent with the 
phenomenological observation that higher-dimensional representations of the flavor group 
are suppressed compared with lower-dimensional representations, and also consistent with 
the additional suppression predicted in the large- A'"c limit of QCD. 
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